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Abstract: We compute the next-to- next-to-leading order QCD corrections to the gluon- 
fusion production of a Higgs boson in models with massive color-octet scalars in the (8, l)o 
representation using an effective-theory approach. We derive a compact analytic expression 
for the relevant Wilson coefficient, and explain an interesting technical aspect of the calcu- 
lation that requires inclusion of the quartic-scalar interactions at next-to-next-to-leading 
order. We perform a renormalization-group analysis of the scalar couplings to derive the 
allowed regions of parameter space, and present phenomenological results for both the 
Tevatron and the LHC. The modifications of the Higgs production cross section are large 
at both colliders, and can increase the Standard Model rate by more than a factor of 
two in allowed regions of parameter space. We estimate that stringent constraints on the 
color-octet scalar parameters can be obtained using the Tevatron exclusion limit on Higgs 
production. 
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1. Introduction 



The Higgs boson is the last undiscovered particle of the Standard Model (SM). The hunt for 
the Higgs boson to uncover its role in electroweak symmetry breaking is being undertaken 
at the Tcvatron, and will soon begin at the Large Hadron Collider (LHC). The CDF and 
DO collaborations at the Tevatron have recently announced a 95% exclusion limit on a 
SM Higgs boson with a mass in the range 162 GeV < rrih < 166 GeV [1], while the LEP 
collaborations have established the limit ruh > 114 GeV [2]. A crucial component of this 
search is the derivation of accurate theoretical predictions for the cross section for Higgs 
production. The dominant hadronic production mechanism, gluon fusion through a top- 
quark loop, is known exactly through next-to-leading order in perturbative QCD [3,4]. In 
the effective theory with rrit — t- oo, both the NLO corrections [5] and the NNLO corrections 
are known [6-8]. When normalized to the full mj-dependent leading-order result, this 
effective theory reproduces the exact NLO result to better than 1% for Mh < 'imt and 
to 10% or better for Higgs boson masses up to 1 TeV [9]. The QCD radiative corrections 
drastically alter the Higgs production cross section prediction; for example, the gluon 
fusion cross section is increased by roughly a factor of three above the LO prediction at 
the Tevatron after the NNLO corrections are included. Only at NNLO is an accurate 
prediction free from debilitating uncertainties obtained. Updated cross sections for Higgs 
production in gluon fusion for use at the Tevatron and LHC are available in Refs. [10,11]. 
For a recent review of the status of theoretical predictions for Higgs boson production in 
the SM, see Ref. [12]. 

The properties of the Higgs boson can be modified in theories with additional particles, 
and measurement of these properties consequently serves as a window into physics beyond 
the SM. In regions of supersymmetric parameter space, the production mechanism gg — >■ 
/i — > 77 can be changed by up to a factor of two [13]. In extra-dimensional theories, mixing 
between the Higgs and the scalar radion can modify the Higgs production rates by orders 
of magnitude [14]. The gluon- fusion mode is particularly sensitive to new states, since 
they can contribute to the cross section at the same one- loop order as the SM particles. 
It is therefore important to accurately predict the gluon-fusion production cross section in 
theories beyond the SM to assure that signatures of new physics are conclusively identified. 

Several interesting extensions of the SM introduce new scalar states transforming in 
the adjoint representation under the color gauge group. It was shown in Ref. [15] that 
a scalar particle transforming as (8, 2)i/2 under SU{3) x SU{2) x U{1) is the only new 
scalar representation with non-trivial electroweak quantum numbers that can couple to 
quarks without introducing additional flavor violation beyond that present in the SM. 
As the scalars can couple to the SM Higgs at tree-level, such particles can modify the 
gluon-fusion cross section through their propagation in loops. The induced shifts to the 
gg ^ h cross section induced by (8,2)i/2 scalars were studied at NLO in Ref. [16], and 
were found to be large. Scalars in the (8, l)o representation can arise in theories with 
universal extra dimensions [17, 18] and in technicolor models [19]. The primary decay 
for such states is expected to be into either bb or tt, depending on the scalar mass and 
other model parameters. The Tevatron experiments can search for these states via pair 
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production of (8, l)o scalars, leading to a four 6-jet final state. The search reach in the 
scalar mass was recently estimated to be 280 GeV [18]. The direct search for these scalars 
is rendered difficult by the large QCD background. Reduction of the background requires 
significant cuts that reduce the signal and therefore the search reach. It is possible that 
indirect searches for these scalars, such as via their influence on the Higgs production cross 
section, can probe masses competitive with direct searches. 

We compute in this manuscript the NNLO corrections to the Higgs boson production 
cross section in models with a (8, l)o scalar. To be as independent as possible from the 
origin of this scalar, we study a simple model that couples this scalar to both QCD and 
the Higgs doublet via renormalizable operators. We utilize the effective-theory approach 
valid when both the SM top quark and the new scalar are heavier than the Higgs boson. 
As explained above, we expect this to be an excellent approximation for Higgs masses in 
the interesting range. We focus here on the (8, l)o state as a first case because it leads to a 
significant technical simplification when handling the effects of the scalar-sector potential. 
As we explain later, the need to include the scalar self-couplings first arises at NNLO. 
We compute the Wilson coefficient describing the interaction of the Higgs boson to gluons 
mediated by both the top quark and the color-octet scalar through NNLO in the QCD 
coupling constant. We derive all renormalization constants required in the scalar sector. 
All results arc presented in compact analytic expressions. We derive the renormalization- 
group equations governing the evolution of the scalar-sector couplings, and use them to 
determine the likely ranges of the various parameters which appear. 

In addition to the analytic results described above, we study the phenomenological 
implications of the color-octet scalar for the Higgs production cross section at both the 
Tevatron and the LHC. Only at NNLO is the scale dependence sufficiently reduced to 
allow precise predictions for the scalar- induced effects to be obtained. The deviations from 
the Standard Model prediction for Higgs production are large at both colliders. Shifts 
in the Higgs cross section larger than the errors coming from scale variations and parton 
distribution functions are obtained for scalar masses approaching 1 TeV. Deviations of 
a factor of two arc obtained for scalar masses near the estimated direct search reach of 
ms ~ 300 GcV. We therefore believe that the current Tevatron exclusion limits on the 
Higgs boson would yield constraints on the color-octet parameter space competitive with 
direct-search constraints. The indirect constraints depend on an undetermined Higgs-scalar 
coupling, and are therefore more model dependent. However, from a low-energy perspective 
there is no symmetry reason to expect a small value for this coupling, and therefore a large 
region of allowed parameter space would be tested. 

Our paper is organized as follows. We describe our model for the (8, l)o scalar in 
Section 2. In Section 3 we describe our calculation, and present analytic results for the 
Wilson coefficient and all required renormalization constants in the scalar sector. Numerical 
results for both the Tevatron and LHC are presented in Section 4. We conclude in Section 5. 

2. Details of the Model 

We begin with the following Lagrangian, which describes the Standard Model coupled to 
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a color-octet scalar in the (8, l)o representation: 

Z:^"" = CsM + Tr [D^SD^'S\ - ml Tr [5^] - gl G^s Tr [S^] ' - XiH^H Tr [S^] 

-A,(^i/t^-^) . (2.1) 

S denotes the matrix- valued scalar field S = S^T^, H indicates the Higgs doublet before 
electroweak symmetry breaking, v is the Higgs vacuum-expectation value, and -D^ is the 
covariant derivative for adjoint fields. The Higgs quartic coupling has been explicitly 
included to define its normalization. It will be used later when deriving the allowed range of 
the scalar couplings. After electroweak symmetry breaking, the Higgs doublet is expanded 
as H = (O, {v + h)/^/2) in the unitary gauge. The mass of the color-octet scalar becomes 
m| = m'g + \\v^ jl. The Feynman rules which describe the scalar couplings to the Higgs 
boson h and to gluons are easily obtained from Eq. (2.1). The free parameters which 
govern the scalar properties are m^, Ai, and G45. We note that higher-order operators 
which break the 5 — ?■ — S* symmetry present above and which allow the scalar to decay are 
obtained in explicit models which contain this state [17,18]. We neglect them here since 
we anticipate that they have little effect on the gg ^ h production cross section. 

Before continuing, we comment on the appearance of the Tr [5^] ^ term in the scalar 
potential. A quartic-scalar coupling is generated by QCD interactions even if it is set to zero 
at tree-level. At NNLO the quartic coupling must be included to obtain a renormalizable 
result due to the contributions of three-loop diagrams such as shown in Fig. 1. We include 
this operator in the tree-level Lagrangian with a coefficient scaled by g^, the QCD coupling 
constant squared, to permit an easier power-counting of loops. 



Figure 1: Example of a three- loop diagram necessitating the inclusion of the quartic-scalar term 
in the Lagrangian. 
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3. Calculational details and analytic results 

We discuss here our derivation of the effective Lagrangian describing the interaction of the 
Higgs boson with gluons through NNLO in the QCD coupling constant. When both the 
adjoint scalar and the top quark are heavier than the Higgs boson, they can be integrated 
out to derive the following effective Lagrangian: 

H 

V 

where Ci is a Wilson coefficient and the operator Oi is 

1 
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C'" = /^"Sd-Ci-0,, (3.1) 



Oi = ^ G'^G'"'"' . (3.2) 
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The effective Lagrangian for the gluons and hght quarks, CJqcd ^ has the same form as 
CJqcd except that the fields and parameters it contains are rescaled by decouphng constants 
that account for the effects of the heavy states. In order to distinguish the fields and 
parameters occurring in the full Lagrangian from those in the effective Lagrangian, we 
denote the latter ones with a prime. In this manuscript we integrate out the top quark and 
the color-octet scalar in a single step. Our calculation is therefore a two-scale problem. 
This is clearly demonstrated by the example diagrams contributing to the NNLO Wilson 
coefficient shown in Fig. 2. If a large hierarchy exists between the scalar and top masses, 
a two-step procedure in which the scalar and top quark are integrated out separately can 
instead be used. 
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Figure 2: Example two-scale diagrams contributing to the Wilson coefficient at NNLO. 



The Wilson coefficient Ci can be constructed by computing the amplitude for gg ^ h 
in the limit that the initial gluon momenta vanish, as reviewed in Ref. [21]. We start from 
the relation 



{N^-l){d-2){p^.p2f 
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0,aia2 



{Pl,P2) 



lpi=P2=0 



(3.3) 



between the bare Wilson coefficient Ci and the bare amplitude for the process gg ^ H in 
the full theory. Here, pi and p2 are the momenta of the two gluons, N is the number of 
colours, and d = 4 — 2e is the dimension of space-time. The factor (^3 is the bare decoupling 
coefficient by which the bare gluon fields in the full and effective theories are related: 
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^3 '^u ' 



(3.4) 



It can be expressed as 



^3 



(3.5) 



where 11^ is the transverse part of the bare gluon self-energy in the full theory. In order 
to obtain the Wilson coefficient through NNLO in the QCD coupling constant, A^Jl'f/i"^ 
is needed up to three loops, while I^Q{p = 0) is needed through two loops. We can set 
all scaleless integrals to zero in dimensional regularization, and as a result, only diagrams 
containing at least one massive scalar propagator contribute to Hq{p = 0) and to the 
right hand side of Eq. (3.3). We perform our calculations in a covariant gauge with gauge 
parameter ^, leading to the following gluon propagator: 



gf^^ + C- 



(3.6) 



- 4 - 



In all diagrams, terms up to first order in ^ are retained. All ^-dependent terms cancel 
in li^ip = 0) up to two-loop order, and in up to three-loop order, demonstrating 

that Ci is gauge independent. 

Initially, all quantities appearing in eq. (3.3) are expressed in terms of the bare masses 
m^, nij,, and the bare coupling constants 5°, G^g. We derive here the renormalization 
constants that relate the bare parameters to physical ones. The scalar and top-quark 
masses are renormalized in the MS scheme according to 

i^s = V^rns 'ms, rriT = Zmj. rnr, (3.7) 

with 



=1-^(9- IOG45) 



+ «^ I 7:7r^T^ — + 



237 - 6 n/ - 360 G4S + 260 Gig -53 + 10 n/ + 480 G45 - 100 G^^ - 72 mf. /m| ~ 



and 



32^2 Me 

(3.8) 



, a 2 /83-4n, -553 + 20n/\ 



We note that the appearance of the m^/rrig term in the expression for Z^s is not un- 
expected. The scalar mass receives an additive contribution from the top quark since no 
symmetry prevents it, and as a result of attempting to force a multiplicative renormal- 
ization in Eq. (3.7) this mass ratio appears in Z^g. Equivalently, one can formulate the 
renormalization by introducing a matrix of constants: 

(m0)2 = ^ Z,,(mO)^ (3.10) 

j=T,S 

with i = T, S. The appearance of the m^/rrig term in Eq. (3.7) is then converted into the 
presence of an off-diagonal term Zst in this matrix. 

The quartic coupling first appears in the gg ^ H amplitude at the two-loop level, and 
its renormalization is therefore only required to one loop: 

G^5 = Z45G45, (3.11) 

with 

We convert the bare strong coupling constant of the full theory into the bare coupling of 
the effective theory using decoupling constants obtained from the ghost self energy and the 
ghost-gluon vertex up to two- loops, as described in detail in [21,22]. We then renormalize 
it in the effective theory using the MS scheme: 

gf = l,'Z'gi, (3.13) 
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where 



ryl -I / ^0 , /2 3 ^0 ^1 

Zi„ = \ — a - — h a -—^ — — 
^ 2e V8£2 4£ 



(3.14) 



Here, 
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(3.15) 



and the first two coefficients of the /3 function are given by 



/3o 
/3i 



1 / 2 
- 11 - -n; 
4 V 3 



16 V 3 ' 



(3.16) 
(3.17) 



The Wilson coefficient itself requires a renormalization factor [23-25]: 



^11 



(3.18) 



with 



= H-a ha — 

Zii £ e 



(3.19) 



We split our final result for the renormalized Wilson coefficient into three contributions: 



Ci = Ctth + CssH + Cts- 



(3.20) 



The known SM Wilson coefficient Ctth arising from just the top-quark contributions is 
given by [22,26,27] 



Ctth = — 



a' 11a' 



12 



12 



+ a' 



/3 



(-2777 + 684 Lt) + — (67 + 64 Lt) ni 

864 ^ 288 ^ 



(3.21) 



The contribution CssH arises from diagrams involving the coupling of the adjoint scalars 
to the Higgs boson. Examples of such contributions are given by Fig. 1 and the left-most 
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diagram in Fig. 2. The full contribution to CgSH is given by 



CsSH = - 



2 I a' 



2m| 1 4 



/2 



33 5 Gis 
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16 8 
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(3.22) 



9L5(-43 + 8x^) 3(76-3895a;^ + 257x'') 
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1024x2 



^ , -705 575 L5 51n(x)\ 3 (76 + 37 + 86 a;^ + 225 x^) 



64 



96 



24 



2048x3 



Li3(x) - Li3(-a;) 



+ ln2(a;) 



-228 + 41 x2 - 192 + 675 3 (76 + 37 + 86 x^ + 225 x^) 

H „ X 



2048 (-1 +x) x2(l + x) 



ln(l + x) — ln(l — x) 



4096x3 



76 - 111x2 + 159x^ 76 + 37x2 + 86x^ + 225x6 , , 

+ 31n(x) <; ^^24^2 ^^KX^?. ( Li2(2;) - Li2(-x, 



2048x3 



The remaining contributions Cts come from the adjoint scalar propagating in loops, but 
where the Higgs couples to the top quark. This piece gives the correction to the Higgs 
production cross section that would result if the scalar-Higgs coupling were set to zero. 
An example contribution is given by the right-most diagram in Fig. 2. We used Ref. [36] 
for the non-trivial 3-loop master integrals with two scales and four propagators. The full 
contribution begins first at the three-loop order, and is given by the following expression: 



,'3 



Cts = a 

+ ln(x) ■ 

+ ln2(x) 



9L5x2 2052 + 1075x2 + 1755 x"" 



8 9216x2 

684 + 409 x2 + 1431 x^ 3 (76 + 37 x2 + 86 x^ + 225 x^) 



(3.23) 
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2048x3 



Li2(x) - Li2(-x; 



228 + 11 X- - 192 + 075x6 3 (7G + 37 .r^ + 8G,r^ + 225x6) 

- H ^ '- X 



2048 (-l+x)x2 (1 + x) 4096x3 

I , , , \\ 3 (76 + 37x2 + 86 x'' + 225x6) 
I ln(l ' ""^ ^^''^ ^^ 1 V. I V 



+ x 



ln(l 



+ 



2048x3 



Li3(x) - Li3(-x) 



In these expressions, we have set Lj = In (mj//i) for i = [S, T), x = mT/ms, and we have 
used niT and ms to denote the MS scalar and top quark masses. 



4. Numerical Results 

We now present numerical results for the Higgs boson production cross section in gluon 
fusion, to study the deviations induced by the color-octet scalar. We include the effects 
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of the top quark, the scalar, and also the bottom quark on the Higgs cross section. We 
comment first on precisely what terms we include in the cross section. The leading-order 
amplitude for the gg ^ h process takes the form 



where the subscripts T, B, S respectively denote the top, bottom, and scalar contributions. 
Upon squaring this amplitude, interferences between the contributions of each particle are 
obtained. We denote by ctt+s the terms obtained by squaring together the top and scalar 
amplitudes, and keeping both the interference term and the pieces from each separate 
particle squared. We let aTB,crsB denote the interferences between the bottom-quark 
amplitude with the top and the scalar pieces, respectively. For the cross section at the n-th 
order in perturbation theory, we use the following expression: 



^EFT denotes the ratio of the n-th order cross section over the LO result, with both 

computed in the effective theory defined in Eq. (3.1). The cross section multiplying the 
i^'-factor is the LO cross section with the exact dependence on the scalar and top-quark 
masses. The remaining terms account for the scalar-bottom interference, the top-bottom 
interference, and the bottom-squared contribution at LO with their exact mass depen- 
dences. In the SM, the scaling of the exact LO cross section by the EFT K-factor is 
known to furnish an approximation accurate to the few-percent level or better for Higgs 
masses below roughly 400 GeV [9]. For the scalar, this approximation has been studied 
against the exact NLO calculation [16], and is again accurate to the 1 — 2% percent level 
for m/i < ms, which is the region we focus on here. If desired, the exact NLO corrections 
to the top-bottom interference and bottom-squared terms [4, 30] can be included, as can 
those to the bottom-scalar interference [16]. These affect the cross section at the 1 — 2% 
level, and for simplicity are neglected. Various cicctrowcak corrections known for the SM 
contribution [10,31-33] are not known for the scalar, and for consistency are neglected. 

We use the the pole mass rriT = 173.1 GeV for the top quark [28], while for the bottom 
quark we use the MS mass with mb(10 GeV) = 3.609 GeV [29]. The choice of pole or MS 
mass for the 6-quark has an insignificant effect on the fractional deviation between the 
scalar-induced cross section and the SM result. We use the MSTW parton distribution 
functions [34] extracted to the appropriate order in perturbation theory. For the scalar 
sector, we must set the parameters Ai, G45, and ms- We perform a renormalization-group 
analysis to constrain the possible values of the scalar couplings. We describe this analysis 
below. 

4.1 Renormalization-group analysis of the scalar couplings 

While the scalar mass ms can be constrained by direct searches performed at the Tcva- 
tron [18], such analyses do not restrict the allowed regions of the scalar couplings G45 and 
Ai. To determine the reasonable range for these couplings, we instead derive constraints 
arising from theoretical consistency of the model. For the potential in Eq. (2.1) to be 
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a" = a^Xsi-mr, ms) K^p^ + '^ssi'ms, mi,) + a^%{mT, mb) + a^%{mb). 



(4.2) 
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bounded at large values of Tr[S'^], we must have G45 > 0. The coupling Ai may be either 
positive or negative. For negative values, the restriction of a bounded potential at large 
field values imposes the constraint |Ai| < 2^/gjG4sXh■ This restriction does not exist for 
Ai > 0. For negative values of Ai the cross section for Higgs production is decreased due 
to destructive interference between the top-quark Wilson coefficient and the scalar con- 
tribution in Eq. (3.22). Since the allowed parameter space for positive Ai is larger, and 
also since Tevatron restrictions on positive Ai are stronger, we focus in this paper on this 
region. 

We derive further constraints by demanding that the couplings do not encounter a 
Landau pole for energy scales up to some cutoff A. At that scale, the model must be 
embedded into a more complete theory that removes the Landau pole. Since the couplings 
increase in the ultraviolet, this restricts their values at lower energies relevant for studies at 
the Tevatron and the LHC. As long as A is smaller than the energies probed by experiment, 
the model defined in Eq. (2.1) provides a consistent framework for the scalar interactions. 

We begin by presenting the system of coupled differential equations which governs the 
one-loop evolution of the MS couplings G45, Ai, and A^. We derive these results from the 
one-particle irreducible Green's functions for the four-scalar interaction, the interaction of 
two scalars with two Higgs bosons and the four-Higgs coupling. The renormalization group 
equations are as follows: 



where L = In/x^, a = g1/{A'K'^), and n/ = 5 is the number of light fermions. The running of 
the strong coupling a including the scalar is well-known. These equations can be checked to 
agree with well-known results available in the literature [20] . For this part of our analysis 
only, we have permitted the Higgs boson to propagate in loops in order to derive the correct 
dependence of these equations on all scalar couplings. 

We solve these equations numerically, starting the evolution at /x = v where the Higgs 
quartic coupling A/i is related to the Higgs boson mass via m\ = 2\h{v)v'^ . We vary the 
starting values of the other couplings \i{v) and Gis{v), and demand that all couplings 
remain perturbative until the scale A = 10 TeV. Stronger restrictions would result if 
perturbativity were imposed up to a higher energy such as the grand-unified scale or Planck 
scale, but we use the weaker constraint here to assure all potentially relevant parameter 
space is included. We show the results of our analysis below in Table 1 for the value 
m/j = 165 GeV, where we fix one coupling to an allowed value and derive the restriction on 
the other. In Fig. 3 we show the maximum allowed value of X\{v) obtained by demanding 
Ai(/Lt) remain perturbative until /x = 10 TeV, as a function of the Higgs mass, for the 
choices G4^s{v) = 1 and G/i,s{v) = 0. 

We note that the combination Aif^/m|. is scale-invariant if only QCD-induced as 
corrections are considered (i.e., the Higgs boson is not allowed to propagate in loops). This 
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dL 
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Fixed Xi{v) 


Fixed Gisiv) 


Xi{v) = 


G^siv) < 1.5 


Xi{v) < 4.3 


GMv) = 


Xi{v) = 2.0 


GAsiv) < 1.3 


Ai(t;) < 3.5 


G^siv) = 0.6 


\i{v) = 4.0 


G4s{v) < 0.3 


Ai(t;) < 2.2 


^45(1^) = 1.2 



Table 1: Restrictions on the couplings Ai and G4S at the input scale fi = v arising from demanding 
perturbativity until /i = 10 TeV. In the three left-most columns, Xi{v) is fixed to a given value, 
and the corresponding constraint on G^siv) is derived. In the three right-most columns, 6*45 (w) is 
fixed to a given value and a restriction on Ai(?;) is found. The value rn/j = 165 GeV is assumed. 
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150 200 250 300 350 400 

rrih (GeV) 

Figure 3: Maximum value of Ai allowed by perturbativity for the choices G4s{v) = 1 and G4s{v) — 
0, as a function of Higgs mass. 

is a consequence of the low-energy theorem for Higgs production [35]. We have verified 
this by explicit calculation in the MS scheme. If instead an on-shell scheme is chosen for 
the scalar mass, one should be able to implement a similar choice for Ai that preserves 
the renormalization-group invariance of the ratio. If that is done, the Wilson coefficient 
presented in Section 3 takes the same form in both the A4S and pole schemes, as the scalar 
mass appearing in and the top mass in Lt can be interpreted as either the MS or pole 
masses. The difference occurs only at NNNLO in perturbation theory, which is beyond the 
scope of our calculation. In presenting our numerical results we interpret the masses of the 
top and scalar as pole masses, which assumes the choice of scheme for Ai described above. 

In our numerical results for the Tevatron, we set Ai to the fixed value Ai = 2.5. For 
LHC numbers we use Ai = 1.5 to allow results for larger Higgs masses consistent with 
perturbativity to be displayed. We assume that the same scheme is chosen for both this 
quantity and the scalar mass so that the ratio is scale invariant. We have checked that the 
contributions from Cts in Ec[. (3.23) are negligibly small, and the effects of the explicit 
mass dependence appearing in L5 are at the percent level. The cross section therefore 
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depends primarily on the scalar parameters through the ratio Xi/rog. Results for other 
values Af, but keeping the mass fixed at m^, can be approximately obtained by studying 
the presented results at a mass m"g given by Ai/m| = Af/ [mPg)"^. We also set G4,s{v) = 1.0, 
and account for the evolution using the leading-order evolution equation in Eq. (4.3). We 
neglect the dependence of G^silA on Ai for consistency with our calculation of the Wilson 
coefficient, where only QCD corrections were included and the Higgs was not allowed to 
propagate in loops. Inclusion of this term has only a small effect on the running of G45. 
The choice of G45 only affects the cross section at the few-percent level. The coupling 
Ai always appears in the ratio Ai?;^/(m|), as is clear from the presentation of the Wilson 
coefficient contribution CssH in Eq. (3.22). Varying the scalar mass therefore adequately 
accounts for possible cross-section deviations obtainable in this model. 

4.2 Results for the Tevatron and the LHC 

We now present numerical results for both the Tevatron and the LHC, for which we assume 
-y/s = 7 TeV. We begin by showing in Fig. 4 the LHC cross section for 777.5 = 300 GeV at 
LO, NLO and NNLO in QCD perturbation theory, to see the effect of including higher- 
order QCD corrections. The renormalization and factorization scales are equated to /ip = 
fiR = ^, and are varied in the range rrih/A < fJ, < rrih, consistent with previous studies 
of the Higgs production cross section in the SM. The scale-variation errors are large at 
both LO and NLO, and the corresponding error bands do not overlap. Only at NNLO can 
a reliable prediction for the cross section be made. It is also clear from this figure that 
large variations from the SM prediction are possible for scalar masses near the expected 
Tevatron limit of ms ~ 300 GeV. The cross section differs from the SM result by more 
than a factor of two for this parameter value. 

To study further the effect of the color-octet scalar on the Higgs cross-section prediction 
at both the Tevatron and the LHC, we show below in Fig. 5 and Fig. 6 the cross sections 
for the example Higgs mass ruh = 165 GeV as functions of the scalar mass at both the 
Tevatron and LHC. Deviations from the SM are visible over scale errors at both colliders 
for scalar masses approaching 1 TeV. The scalar contributions to the Higgs production 
cross section are large, and direct searches arc hindered by the need to pair produce the 
scalars and by the large QCD background. This suggests that the indirect constraints on 
the scalar parameter space coming from the Tevatron Higgs exclusion limit could be as 
strong as the direct search reach. 

It is interesting to speculate on the potential of the Tevatron experiments to exclude 
color-octet scalar parameter space using the established exclusion limit on Higgs boson 
production. As the Tevatron limits include contributions from Higgs production in as- 
sociation with gauge bosons, and also from weak-boson fusion production of the Higgs, 
it is not possible to simply scale current exclusion limits by the ratio of the gluon-fusion 
production cross section with the scalar over that of the SM. The Tevatron collaborations 
plan to release limits utilizing only the gluon-fusion mode to make this scaling possible [37]. 
However, it is likely that scalar masses in excess of the expected direct-search limits can 
be probed. As discussed in the previous section, the scalar contributions depend primarily 
on the ratio Aif^/(m|). It would be interesting to derive the allowed region in the plane 
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Figure 4: Higgs production cross section in gluon-fusion at the LHC for ms = 300 GeV as a 
function of the Higgs boson mass. The bands indicate the scale variation mft,/4 < fi < rah. From 
bottom to top, the bands indicate the variations of the LO, NLO, and NNLO cross sections. All 
other parameters are as described in the text. For orientation, the SM result at NNLO for the 
central value = mhjl is shown. 



of this ratio versus the Higgs boson mass. One caveat is that the Higgs production cross 
section depends on Ai/m|, while the pair production mode used in the direct search only 
depends on ms- The indirect constraint discussed here is therefore more model dependent. 
Nevertheless, since there is no symmetry argument which suggests that Ai is small, the 
Higgs exclusion limit probes a large and relevant region of scalar parameter space. 

To facilitate studies of color-octet scalar effects on the Higgs boson production cross 
section, we display in Table 2 and Table 3 the cross sections for Higgs production for the 
example points ras = 300 GeV and ms = 800 GeV at the Tevatron and the LHC. We 
include in this table the errors arising from both scale variation, and also from imprecise 
knowledge of parton distribution functions (PDFs). For the PDF errors, we show the 
90% confidence-level error estimates derived using the MSTW error eigenvectors. Other 
possible choices for the PDF uncertainty estimate could be used: either the 68% confidence- 
level error, or the combined PDF-|-as error [38]. In each table we also show the fractional 
deviation of the cross section from the SM prediction, which we denote as 5. The deviations 
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rrig (GeV) 

Figure 5: Higgs production cross section in gluon-fusion at tire Tevatron for nih = 165 GeV as 
a function of tlie scalar mass. The bands indicate the scale variation rah/ A < fi < rrih. All other 
parameters are as described in the text. Also shown is the SM cross section with its corresponding 
scale uncertainty. 

are roughly a factor of two for rag = 300 GeV and 20% for ms = 800 GeV. At the Tevatron 
both the scale and PDF errors are roughly ±10%. At the LHC the PDF error reduces to 
a few percent. For both mass points the scalar-induced deviations from the SM are larger 
than the combined errors. 

5. Conclusions 

In this paper we have derived the NNLO Wilson coefficient which describes the effect 
of an (8, l)o scalar on the Higgs-gluon effective Lagrangian. We have presented simple 
analytic formulae for the Wilson coefficient and all required renormalization constants for 
use in future studies. Our calculation revealed an interesting technical feature of scalar 
contributions to the gg ^ h cross section which first appears at NNLO. At this order, 
the quartic-scalar potential must be included in order to properly renormalize the Wilson 
coefficient. In our example with a single (8, l)o scalar, only one additional Lagrangian 
term was required. In theories with a more complicated scalar potential, such as that for 
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ms (GeV) 

Figure 6: Higgs production cross section in gluon-fusion at the LHC for nih = 165 GeV as a 
function of the scalar mass. The bands indicate the scale variation mh/A < /i < rrih. All other 
parameters are as described in the text. Also shown is the SM cross section with its corresponding 
scale uncertainty. 



a (8,2)i/2 scalar [15], the addition of multiple operators is in principle required to obtain 
a finite NNLO result. 

We have performed a renormalization-group analysis of the scalar-sector couplings to 
determine the theoretically consistent region in which perturbation theory can be applied. 
We have also studied the phenomenological predictions for both the Tevatron and LHC. The 
scalar-induced effects are large at both colliders, and are visible over both scale and PDF 
errors for scalar masses approaching 1 TeV. This suggests that the scalar parameters can 
be stringently constrained using the exclusion limit on the Higgs boson production cross 
section obtained by the Tevatron, and that the limits would be competitive with those 
obtained from direct searches for scalar pair production. If the gluon-fusion contribution 
to the Tevatron exclusion limit were separately presented, it would be possible to rescale 
these results to obtain bounds on the Ai/m| combination of scalar parameters. 

In summary, color-octet scalars are an interesting and phenomenologically-rich example 
of physics beyond the SM. The NNLO calculation presented here quantifies precisely the 
effect of a (8, l)o scalar on the Higgs production cross section in gluon fusion, and our 
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Table 2: Table of cross sections for ms = 300 and 800 GeV at the Tevatron. Values are shown for 
Higgs masses between 115 and 200 GeV, and for Ai = 2.5.. The scale and PDF errors are denoted 
by the labels sc and pdf, respectively. The fractional deviations of each cross section from the SM 
for the scale choice fi = are indicated by 5. 

phenomenological study shows that the deviations induced by the scalar are large at the 
Tevatron and the LHC. The exclusion limit on Higgs production set by the Tevatron 
collaborations allows for the scalar parameters to be stringently constrained, and we hope 
that this analysis is undertaken. 
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'^•^^-0.23(sc)-0.14(pd/) 


n 1 nn 
U.iUU 


oOO 


. ,^+0.39(sc)+0.24(p4/) 
^•'J^_0.47(sc)-0.26(pd/) 


n on 

u.yu 


Qp;p; 
oOO 


„ „„+0.18(sc)+0.14(pd/) 
^•"^-0.24(sc)-0.15(pd/) 


n no/1 

u.uy4 


ODO 


. c„+0.39(sc)+0.25(pd/) 
*-^^-0.47(sc)-0.27(p4f) 


n 87 

U.o / 


000 


„ „ ,+0.18(sc)+0.14(pd/) 
Z.D4_Q 24(sc)-0.16(p(i/) 


u.uy i 


375 


^ „q+0.37(sc)+0.24(pd/) 
^•'3y„o.45(sc)-0.27(p4f) 


0.85 


375 


„ ^/,+0.17(sc)+0.14(prf/) 
z.oo_Q 24(sc)-0.16(prf/) 


0.088 


385 


, , q+0.36(sc)+0.24(pd/) 
*-^^-0.43(sc)-0.26(p4f) 


0.84 


385 


„ ,„+0.16(sc)+0.14(pd/) 
^•^'-0.23(sc)-0.15(prf/) 


0.086 


395 


o Q^+0.34(sc)+0.23(p(#) 
'^•^"-0.41(sc)-0.25(p(i/) 


0.83 


395 


„ „ ,+0.15(sc)+0.14(prf/) 
^•'^^-0.22(6-c)-0.15(p4f) 


0.085 


405 


o y-| +0.31(sc)+0.22(pd/) 
''•'-'■-0.38(sc)-0.24(pd/) 


0.83 


405 


„ -,Q+0.14(6-c)+0.13(p4f) 
^•-'-^-0.20(6-c)-0.14(pfl/) 


0.084 


415 


o ,c;+0-29(sc)+0.21(pfi(/) 
^•^^-0.36(sc)-0.23(pd/) 


0.84 


415 


„ P,^+0.13(sc)+0.12(pc(/) 
^•^^-0.19(sc)-0.13(p(i/) 


0.082 



Table 3: Table of cross sections for ms = 300 and 800 GeV at the LHC. Values are shown for 
Higgs masses between 115 and 200 GeV, and for Ai = 1.5. The scale and PDF errors are denoted 
by the labels sc and pdf, respectively. The fractional deviations of each cross section from the SM 
for the scale choice = are indicated by 5. 
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